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18 The accelerating loss of mass from the Greenland ice sheet is a major contribution to current 

19 sea level rise 1 . Increased meltwater runoff is responsible for half of Greenland’s mass loss 

20 increase . Surface melt has been increasing in extent and intensity, setting a record for surface 

21 area melt and runoff in 2012 . The mechanisms and timescales involved in allowing surface 

22 meltwater to reach the ocean where it can contribute to sea level rise are poorly understood. The 

23 potential capacity to store this water in liquid or frozen form in the fim (multi-year snow layer) is 

24 significant 4 , and could delay its sea-level contribution. Here we describe direct observation of 

25 water within a perennial firn aquifer persisting throughout the winter in the southern ice sheet, 

26 where snow accumulation and melt rates are high. This represents a previously unknown storage 

27 mode for water within the ice sheet. Ice cores, ground/airborne radar and a regional climate 

28 model are used to estimate aquifer area (70 ± 10 x 10 km ) and water table depth (5-50 m). The 
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perennial firn aquifer represents a new glacier facies to be considered in future ice sheet mass 
and energy budget calculations. 

The mass of liquid or refrozen meltwater that could be stored in fim pore space throughout 
the percolation zone of the entire ice sheet is estimated to be between 322 and 1,289 Gt 4 . In the 
western part of the ice sheet, the possibility of liquid water persisting within the upper ~10 m of 
the snow/ fim 5 or in moulins 6 during winter is suspected. Discharge measurements at ice marginal 
streams indicate winter water release, suggesting that some meltwater may be stored englacially 
or at the bed and is released months after the end of the melt season 7 . However, there has been 
no account of directly observed liquid water in the firn persisting through the winter on the 
Greenland ice sheet. 

In April 2011, prior to seasonal surface melt onset, the Arctic Circle Traverse (ACT) 
expedition drilled into a liquid water layer in the upper 10 to 25 m of the firn in southeast 
Greenland. The ACT field party extracted four firn cores at sequentially lower elevations on the 
southeast coast (Fig. 1), where in situ snow accumulation observations were previously 
nonexistent. Below 1600 m in this area, spatially and temporally averaged accumulation rates of 
1-4 m w.e. a' 1 are simulated by observationally-constrained regional climate models 8,9, 10 . On 30 
April 201 1, at ACT1 1-A2 (1559 m a.s.l.), a 10 cm diameter firn/ice coring drill extracted a ~1 m 
core segment from 10 m depth that was saturated with liquid water (Fig. 1). The following day, 3 
km to the east at ACT1 1-A (1589 m a.s.l.), liquid water was found at 25 m depth using the drill. 
The thickness of the water layer could not be measured because the drill is not designed to 
operate in water. Air temperatures were -15 °C during drilling. During spring 2011, temperatures 
were below average and surface melt in the area did not commence until June that year 1 ’. 
Therefore, the liquid water found in the fim no doubt persisted throughout the winter. The other 
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two ACT cores were extracted at higher elevations nearby (1806 and 2081 m a.s.l.) and revealed 
no liquid water to the full depth of the 61 m drilling (Fig. 1). 

Ground penetrating radar (GPR, Supplementary Infonnation) profiles were completed 
between the core sites, as well as 10 km below the lowest site. A strong contiguous return 
horizon persists over the lower 25 km portion of the transect (Fig 2). The horizon undulates 
between depths of 9 - 25 m and matches the depth of the water layer top found at both core sites 
to within < 1 m (the precision of identifying the depth is limited to the 1 m length core sections 
drilled). We are thus confident that the GPR is tracing the top of the water layer. 

The top of the water layer cuts across intermediate GPR horizons (Fig. 2), usually 
interpreted as corresponding with annual or event accumulation layers ’ . Below the water 
layer horizon there are no coherent GPR horizons, which can be expected as minimal energy is 
returned from below a strong reflector with a high permittivity contrast such as water 14 . The 
bright horizon gradually fades at the 25 km location at a depth of ~27 m (Fig. 2), revealing 
internal firn layers to depths of —50 m that are traceable up-glacier over the next 82 km to cores 
ACT1 1-B and C (Fig. 1). 

The NASA Operation Ice Bridge (OIB) airborne accumulation radar (AR) 15 overflew the 
core sites and the ground traverse GPR transect 1 1 days prior to the core drilling. A strong 
reflecting horizon is evident at the same location in the GPR transect [Supplementary Fig. SI]. 
The depth to the bright horizon from the GPR and AR agree within 2 m over the 25 km transect 
and the undulations are very similar [r = 0.95, Supplementary Fig. S2]. The depth differences 
may be attributed to lateral discrepancies in the transect locations (< 200 m), and differences in 
radar foot print size and radar frequency. Similar to the GPR data, the AR returns no obvious 
internal layering below the bright horizon [Supplementary Fig. SI]. Based on the depth 
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agreement and high correlation we conclude that the AR is capable of mapping the presence and 
depth to the top of the water layer within the firn. 

In 201 1, prior to melt onset, NASA OIB AR gathered 40,5 12 km of horizontal flight line 
data over the GrIS (March 29 — May 16). All of these flight lines were examined for the 
presence of the water layer. It was identified and manually digitized in 843 km of these flight 
lines acquired between April 8-26, 2011 (Fig. 1). The water layer locations are concentrated in 
the southeast, but are evident in isolated locations in the south and southwest and on the Geikie 
Plateau (near 70° N, 25° W). The mean depth of the water layer top is 23 m with a range of 5 to 
50 m [Supplementary Figs. S3 and S6]. In general the depths are smaller in the southwest 
compared to the southeast, but are influenced by local surface slope (Fig. 2), similar to terrestrial 
groundwater and fim aquifers on temperate glaciers 16 ' 17 . Thus, we refer to this liquid water 
reservoir that persists throughout the winter as a perennial firn aquifer (PFA). Since the radar 
signal is not returning from below the top of the PFA, there is currently no direct measure of its 
thickness. 

The spatial distribution of the PFA suggests that its formation is associated with areas of 
sufficient surface melt coupled with high accumulation. We compare the PFA extent as 
determined by AR results with gridded climate fields. Here we use the accumulation field from 
the Calibrated Polar MM5 10 (Fig. 3) and output fields of accumulation, melt and rain from the 
regional atmospheric climate model RACM02 [Supplementary Fig. S5]. RACM02 includes an 
interactive (with the atmosphere) snow/firn/ice model . In areas where the PFA is found, the 
mean accumulation rate is 1.24/2.22 m w. e. a' 1 (10 and RACM02, respectively). However, there 
is significant variability in the range of accumulation rate associated with the PFA 
[Supplementary Fig. S6]. (Supplementary Information for discussion of the differences in these 
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accumulation grids) Areas of high accumulation are found predominately in the southeast, but 
three other areas in the south and southwest are also identified with local accumulation maxima 
and also contain a PFA (Fig. 3). 

The spatial pattern of the PFA extent for April 2011 simulated by RACM02 is very similar 
to the AR results, with a concentration in the southeast and the three areas in the south and 
southwest (Fig. 4). In RACM02, a combination of high accumulation (> 800 mm yr' 1 ) and a 
large liquid water production (snow melt plus rain > 650 mm yr' 1 ) are necessary conditions for 
PFA formation [Supplementary Fig. S7]. Because RACM02 performs well in detecting the PFA, 
we use it as a first order estimate of PFA spatial extent. Since the model lacks treatment of 
potentially significant firn processes, i.e. inhomogeneous vertical water flow (piping) which 
moves water to depth through cold snow/fim 5 ' 18, 19 and horizontal water flow, which as Fig. 2 
shows is an important process, an estimate of PFA volume is not given here. For uncertainty in 
the modeled PFA extent we use the variability in annual minimum extents during 1992-2011, 
which are typically reached in late April at 70 ± 10 x 10 km , [Supplementary Fig. S8]. After the 
onset of the melt season, extent of liquid water in the Greenland firn sharply increases, to reach 
500 ± 250 x 10 km in July, after which a gradual decrease is simulated. 

The formation process of the PFA is not completely understood; however, its spatial 
correspondence with high accumulation and melt rates (Fig. 3, Supplementary Fig. S7) leads to a 
general hypothesis intended to explain the broad pattern of PFA location. From the RACM02 
model results, we propose that high accumulation insulates the melt season’s liquid water layer 
within the firn from the cold season atmosphere, thereby preventing complete refreezing. This 
allows liquid water to persist throughout the winter until the next melt season, when the PFA 
may be recharged [movie SI]. 


121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 


6 


1 o 

The PFA represents a new glacier facies and a previously unidentified liquid water 
reservoir. Its location in the southeastern ice sheet is consistent with few surface lakes, compared 
to other sectors of the ice sheet . The narrow (< 30 km) ablation area with minimal bare ice area 
in the southeast, due to the relatively high accumulation gradients and therefore steep ice slopes, 
does not accommodate lake basin formation with accompanying supra glacial stream networks 
and moulins delivering water toward the bed as prevalent along the western portion of the ice 
sheet. In its place, the deep fim layer provides an alternative liquid water reservoir in winter, 
which may exceed the mass of liquid water stored in supraglacial lakes. This contrast in liquid 
water storage mechanisms implies that surface mass balance, thermal properties of the ice, and 
effective water pressures at the bed and consequently ice dynamics in the southeast are likely 
very different from those in the more extensively studied western and northern Greenland ice 
sheet. The persistence of liquid water in the fim also has implications for the ice sheet energy. If 
atmospheric warming ceases, refreezing the liquid water requires a significant amount of 
additional energy before the firn layer can start to cool. The PFA could thus represent an 
increasingly important mass and energy reservoir, as both melt and accumulation on the 

Greenland ice sheet have increased in the past and are projected to increase in a future warming 

22 

climate" . 

Methods 

GPR data were collected with a Geophysical Survey Systems, Inc. SIR-3000 controller and 
a 400 MFIz center frequency antenna. The vertical resolution in firn is 35 cm ’ , finer than the 
annual layering in this area. The sampling was set at four traces per second, 2048 samples per 
trace. To increase the signal-to-noise ratio an initial stacking of six traces was performed. Post 
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processing horizontal spatial smoothing involved averaging an additional eight traces, to increase 
the signal-to-noise ratio and minimize the influence of cm-scale vertical ice pipes or channels 
present in the percolation zone . A time dependent gain was used to compensate for signal 
attenuation within the fim. The maximum two-way-travel time (TWT) range was set to 500 ns, 
yielding a -0.24 ns sample interval and allowing scanning of the top -46 m of the fim. Because 
the top of the water layer is mostly found in the upper fim column, radar TWT is converted to 
depth below the surface assuming a constant electromagnetic wave travel at 1.94x10 m s' into 
the fim. This travel velocity corresponds to a depth-averaged firn density of 650 kg m' . We 
compare this method with the TWT-depth conversion described by (13) using the relationship 
between velocity in the fim and ACT1 1-A firn density profile 24 . The difference between the two 
methods does not exceed 50 cm for the first 20 m of the firn column, and with the lack of 
detailed density profiles (despite at our fim-core locations), we favored the first method in our 
analysis. 

The GPR did not have an integrated GPS, therefore a roving GPS unit was attached to the 
snowmobile towing the GPR sled and collected a point every five seconds.. GPS data processing 
was done using the on-line Canadian Spatial Reference Service - Precise Point Positioning. This 
processor uses GPS orbit and clock information to enhanced positioning precisions in the 
International Terrestrial Reference Frame via a kinematic processing mode. To geo-reference the 
final GPR radar data, the processed GPS data were used by matching the GPS time to the starting 
point of each GPR radar image in post processing to yield a 10 cm-scale topographic profile 
coincident with all GPR lines. A linear interpolation of the 5 s GPS points was made to obtain a 


GPS coordinates for each GPR trace. 
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The Accumulation Radar 15 (AR) is a combined stepped-chirped system built by the Center 
for Remote Sensing of the Ice Sheets (CReSIS), operates from 550 to 900 MHz when flown on a 
P3 aircraft typically 500 m above surface with a vertical resolution in ice of 28 cm over an 
effective footprint of approximately 30 m (https://www.cresis.ku.edu/). All of the AR radar 
images (example: fig. SI) from the NASA Operation IceBridge (OIB) flightlines were manually 
inspected for presence of a water layer representing the top of the PFA. This was characterized 
by a strong subsurface horizon with no internal layers below. The top of the PFA along with the 
snow surface was screen-digitized on the corresponding radar images. The time difference 
between the surface and reflection horizon was converted to depth to the top of the PFA using 
the same wave velocity of 1.94x10 m s' as was used for the GPR depth calculations. A direct 
comparison between the depth to top of PFA derived from the GPR and AR is made for the 25 
km segment centered on the cores that drilled to water (ACT-1 1A and ACT-1 1A2, fig. S2). The 
depth to top of the PFA along the OIB flight lines is shown in fig. S3. 

While the spatial patterns between the PFA as mapped by the AR and simulated by 
RACM02 are similar (Fig. 4), potential reasons for their differences (apart from the obvious 
uncertainty in the model results) are described here. The areas of RACM02 simulated PFA that 
are not mapped as PFA in AR flight lines may be due to several reasons: 1) water is present but 
subsurface returns from the radar are not detectable because strong surface returns (clutter) from 
rough crevassed surfaces mask the weaker water layer return, 2) water is present but subsurface 
returns from the radar are not detectable due to subsurface clutter (refrozen ice bodies above the 
water layer), 3) liquid water is not present due to internal drainage through crevasses, a process 
not included in RACM02. Locations where surface crevassing is observed coincident with 
simulated RACM02 PFA that is not detected by AR are shown in Fig. S4. These three scenarios 
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189 could explain the lack of mapped water layer along flight lines with modeled water over the 

190 lower elevation portions of the numerous outlet glaciers along the south east coast. The contrary 

191 situation, with mapped PFA that is not simulated, is isolated to areas in south west and could be 

192 associated with model uncertainty and limited model resolution. 
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Fig. 1. Perennial Firn Aquifer locations on the Greenland ice sheet detected by radar and firn 
cores in April 2011. NASA’s Operation IceBridge airborne Accumulation Radar flight lines are 
gray and locations of detected PFA are magenta dots. The red line represents the Arctic Circle 
Traverse 2011 with PFA fim-core locations and names (blue diamonds) and dry firn core 
locations red diamonds). The green line corresponds to the Arctic Circle Traverse 2010 that 
found no PFA evidences from firn cores (green diamonds). The ice sheet margin is blue and the 
black segment on ACT-1 1 line (inset) matches the GPR echogram (Fig. 2). 


Fig. 2. Profile of the top of the PFA from ground penetrating radar along ACT-1 1 traverse 
including PFA firn-core locations (ACT1 1-A and ACT1 1-A2). a, Surface elevation profile from 
simultaneously acquired GPS and topographically corrected GPR PFA top horizon. This 
indicates the depth to top of the firn aquifer is influenced by the local topographic slope, b, GPR 
echogram with the top of the firn aquifer as the bright contiguous horizon cutting the numerous 
internal firn reflecting horizons. Location of the GPR profile is shown in Fig. 1 . 


Fig. 3. Annual snow accumulation (1958-2008) from regional climate model with output 
calibrated by ice core values (color). Terrain elevation contours are white. NASA Operation 
IceBridge flight lines are gray. The ACT-1 1 traverse is red. Locations of radar-retrieved firn 
aquifer positions from the OIB Accumulation Radar are illustrated as black dots. 


Fig. 4. Modeled liquid water content (LWC) in the firn and detected PFA from airborne radar. 
The simulation of LWC is from RACM02/GR for April 2011 (color). OIB flight lines (gray), 
ACT-1 1 traverse (red) and locations of PFA from OIB radar (black dots) are all data acquired in 
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297 April 2011. The LWC is integrated for the entire fim column from the surface down to 

298 approximately 20 m, varying with location (see methods for details). 
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